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Abstract 

In a frame of the t — J — G model we derive the microscopical expression 
for the circulating orbital currents in layered cuprates using the anomalous 
correlation functions. In agreement with /z-on spin relaxation (//SR), nuclear 
quadrupolar resonance (NQR) and inelastic neutron scattering(INS) experi- 
ments in YBa2Cu306+x we successfully explain the order of magnitude and 
the monotonous increase of the internal magnetic fields resulting from these 
currents upon cooling. However, the jump in the intensity of the magnetic 
fields at T c reported recently seems to indicate a non-mean-field feature in 
the coexistence of current and superconducting states and the deviation of the 
extended charge density wave vector instability from its commensurate value 

Q~(7r,7r) in accordance with the reported topology of the Fermi surface. 
74.72.-h, 74.25.-q, 74.20.Mn, 74.25.Ha 
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A possibility for a staggered orbital current phase formation in layered cuprates has 
attracted much interest recentlyllhS. Most importantly, it was shown that most of the 
observed properties referred to a pseudogap phenomenon can be naturally explained in an 
extended charge density wave (CDW) scenario with a complex order parameter phase forma- 
tion (shortly s + icf-CDW) in underdoped cuprates. The real s— wave symmetry component 
corresponds to a conventional charge (or spin) density waves whereas the imaginary part 
of the order parameter has a d— wave symmetry and corresponds to the staggered current 
phase. Different kind of experiments can be interpreted in favor of the staggered orbital 
current phase such as an observation of the orbital antiferromagnetism in YBa2Cus06+ y by 
means of inelastic neutron scattering(INS) experiments reported in Refs.!'! and zero-field 
muon spin relaxation (/xsR) experiment si. Moreover, recent investigations using nuclear 
magnetic resonance (NMR) technique indicate the presence of the internal fluctuating mag- 
netic fields in the superconducting state of layered cuprates§~0. Most importantly, the 
observed enhancement of the magnetic moment's intensity at T C B! seems to indicate an in- 
trinsic and a non-trivial relation between the superconducting and the pseudogap phases. In 
this connection a microscopical analysis of the network patterns and the strength of orbital 
currents becomes very actual. 

In general, the possibility of the s + id-CDW phase formation is related to a divergence of 
the dynamical charge susceptibility at wave vector Qj « (tc, tc) in the first Brillouin zone and 
was demonstrated recently for cuprates@. Here we derive the analytical expression for the 
current flow and show how its orbital contour can be reconstructed for any arbitrary chosen 
instability wave vector Qj. Most importantly, we calculate the intensity of the resulting 
internal magnetic fields and the corresponding orbital magnetic moments. We find that 
its enhancement at T c may result from the presence of a relatively small s— component 
of the extended CDW. The latter agrees well with an observation of the increase of the 
NQR linewidth at Cu(2) site (see Ref.0). In addition, the non-mean-field character of the 
coexistence of superconductivity and s + id-CD W phases has to be taken into account. 

Hamiltonian and general expression for the current flow. In our analysis we start from 
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the following t — J — G model Hamiltonian: 



H = Y,u^ d ^r d + T,-h (SiS,)-^ +E^M- , (i) 



i>3 



d_ 

dt 



where =\ i,a >< i,/3 | are projecting Hubbard-like operators. Symbol pd corresponds 
to a Zhang-Rice singlet formation with one hole placed on the copper site whereas the 
second hole is distributed on the neighboring oxygen sitesEl. Here £y is a hopping integral, 
Jij is a superexchange coupling parameter of copper spins, and cr = ±1/2. 8^ = \I/f ' p is 
a hole doping operator. As in Ref.@ we also use the parameters of a screened Coulomb 
repulsion of the doped holes at different sites, Gij. The quasiparticle energy dispersion and 
the correlation functions were calculated in a Roth-type of a decoupling scheme for the 
Green's functions@'0. 

The network patterns and the strength of the orbital currents can be obtained using the 
charge conservation low 

J ' pdV = J jdS . (2) 
The operator of the fluctuating charge per unit cell with number i is given by 

ctPfH* = c* = $ , (3) 
that obeys the equation of motion 

ih^ = [SiH] . (4) 

Calculating the commutator with Hamiltonian ([I]) we arrive to the following expression 

ih^St = e £ U^f>°%* d - e £ t^f> ^ , (5) 

where the right-hand side of this equation is a current operator. In order to calculate its 
thermodynamic value along the link < ij > we make the Fourier transform of Eq. fl5|). Then 
the probability of the hopping from i to j site can be written as: 

< yf^'P* >=1^< m^m^f d > expHkR, + zk'R,-) , (6) 
N k,k' 



whereas of the reverse process is given by 



N 



(7) 



k,k' 



Since the hopping integral is a real quantity the current flow will be proportional to the 
difference of Eqs. (|5|) and (|7|): 



< \^P d ' a ^' pd — \$j pd > a \$ CT > pd >- 

1 3 3 1 



< K d ' a y a v d > {exp(-ikRi + ik'Rj) - exp(-ikR j + zk'R,)}. . 



N 



(8) 



k,k' 



At T < T* one have the following non-zero expectation values: < ^/ Pd ' cr ^)'^ pd > ? < 
*fiq^k Pd > and < *k^k+Q >• Since the first one does not contribute, we have 

< \$ pd ' a \$ a ' pd — \$f pd > a \$f a ' pd >— 

% 3 3 1 

T7 E < *k+q*k Pd > {exp[-i(k + Q)Ri + ikRj] - exp[-z(k + Q)R, + zkR;]} + 
^ k,k' 

+T7 E < ^X'+q > {exp[-zkR, + z(k + Q)Rj-] - expf-zkRj- + z(k + Q)R,;}. (9) 
^ v k 

In our case the pseudogap order parameter is expected to be the complex (s + z<f)-CDW. 
Therefore, it is useful to separate the correlation functions into two parts: Re < ^k+Q^k^ > 
and Im < ^k-fQ^k^ > - ^ is straightforward to right further as 

» j j i 

2i 



— E Im < K%^ P > {cos[kR, - (k + Q)Ri] - cos[kR - (k + Q)R,-]} 



2i 



+9 E ^ e < *k+Q*k Pd > {sin[kR, - (k + Q)Rj] - sin[kR, - (k + Q)R,-]}. 



N 



(10) 



For the lattice with a mirror plane symmetry perpendicular to the x and y axis the integrals 
over the first Brillouin zone containing sinkRjj vanish. Thus, in the functions cos[kRj — 
(k+ Q)Ri]-cos[-kR + (k + Q)R i ], and sin[kR j -(k + Q)R. ; ]-sin[-kR + (k + Q)R j ] one 



coskRjj and 



sin QR, — sin QR ; 



coskR 



3H 



can leave only their parts cos QRj — cos QRj 
respectively. Then, the contribution to the current flow along the x- axis due to the nearest 
hopping can be calculated: 



J' 11 = [cosQRj - cosQRj] Y. Im < *Sq*k'"' > coskR 



sin QR^ - sin QRJ Y, Re < ^Sq^k^ > coskRy. 



(11) 



Let us discuss at the beginning the simplest case of Q = (ir, ir). One can immediately see that 



the second term of Eq. ( |TT| ) vanishes. Eq. ( |TT| ) allows easy to display the network patterns 
for the different symmetries of the order parameter (s-, d- and so on). Most importantly, 
for the pure d-wave symmetry order parameter Im < ^k+Q^k^ >~ cosk x — cosk y (see 
also Ref.0) and hence the current network pattern is directly mapped on the well-known 
flux-phase state0. 

In general case the period of the current pattern is given by a x = 2tc/Q x , a y = 2n/Q y . 
Note, there are also other contributions to the network patterns due to the next-nearest- 
(t 2 ) and next-next-nearest- (£3) neighbors hopping. These parameters are needed for the 
describing the real Fermi surface (see for example Ref. [18|). The contribution due to t 2 is 
given by 



J( 2 ) = U 2 V2^- [cos QR, - cosQRj £ Im < #£ Q #£ pd > coskR^+ 



h n 

-2 

: iv 



+^t 2 V2^ T [sin QR^ - sin QR^ £ Re < > coskR 



(12) 



and due to U as 



j( 3 ) = ±t 3 — [cosQR; - cosQR,-] £lm < K%^ pd > coskRy+ 



sin QR^ - sin QRj ^ Re < ^f Q ^ pd > coskRy 



(13) 



Note in Eq. (|I~2|) indexes i and j refer to the next-nearest neighbors whereas in Eq. ( |13"D i 
and j refer to the next-next-nearest neighbors. 

The required correlation function can be calculated straightforwardly in a mean-field- 
approximation and is given byEl: 



j P d,a j <T,pd PpdPk 



E 



tanh 



ifc 



E 



Ik 



+ 



2k B TJ E, 



tanh 



-2k 



E 



2k 



2k B T 



+ 



pd 



X 



G 

— (e k -e k+Q ) + A k A* k (G k + G* k ) 



1 



5 



■ tanh 



( E 



lk 



\2k B T) R 



1 



tanh 



2(R lk - E 2k y 

E 2 k 



x 



>2k 



2k B T 



(14) 



where 



E 2 lkak = (£fc+ ! fc+ ° )2 + A fe A* + G k G* k ± 



±- 

2 



[4 - e 2 k+Q f + 4G*G* fc ( £fc + e k+Q f + A fe A*(G fe + G* k ) 2 ^ . (15) 



Here G k = S{T) + tG {T)(cosk %CL cos ky& ) is an extended CDW gap and = 
° 2 (cos fc x a — cosfc y a) is a superconducting <i-wave gap. It is important to note that 
the real part of the correlation function (i"3) contains the term which is proportional to the 
superconducting gap AfcA£. Therefore, at T=T C the value of the correlation function should 
display a 'step' which can be responsible for the corresponding 'jump' of the current flow. 
Keeping in mind this qualitative idea let us now turn to the numerical calculations. 

For Q = (tt, 7r) the resulting current is manly determined by Eq. fllTD . Using the equation 
for the id — CDW order parameters in a mean field approximation (see for details Ref.0) it 
is straightforward to prove that 

E /m < Vg&K** > ^sk x a « + S). (16) 

k J\ + ZCri 

Here J\ t and G\ are parameters of the superexchange and the screened Coulomb repulsion 
of the holes at the nearest copper sites taken to be 120 meV and 135 meV respectively. 5 
is number of holes per one unit cell. Comparing Eq. ( |T6| ) and Eq. flTT|) one sees that the 



temperature dependence of the current strength is almost the same as for the order 
parameter Gq(T). The latter was calculated self-consistently in Ref.El and as was shown the 
results are sensitive to the details of competition between superconducting state (SC) and 
id-CDW state. Below T c the superconductivity tries to push out id-CDW state and as a 
consequence the order parameter G (T) goes down at T < T c . We shall describe this effect 
approximately as 



G (T) = Ak B T c x^l - (T/T*) 2 - y6(T c - T)k B T c ^l - (T/T c ) 2 (17) 
where x and y are the parameters of coexistence of the superconducting and the s + id-CDW 



states. Here A (T) = 2.4/c B T c i/l — (T/T c ) 2 , 6{T C — T) is a usual theta function and x is 



a parameter that depends on a doping level. According to the mean field calculations^! at 
T* « 3T C x is equal one and near the optimal doping i.e when T* « T c , a; ~ 1/4. 

In more general case one expects that (7r,7r). Then the total current is a sum of 
+ J^ 2 -* + J^ 3 -* and the relation between the current flow and the gap is not as simple as 
in Eq.flT7]). The deviation of Q from (tt, 7r) is naturally expected as a consequence of the 
changes in the topology of Fermi surface away from half-filling. This is also can be seen 
from our previous analysis of the dynamical charge susceptibility which becomes divergent 
along the contour around Q = (ir, ir) (see for details Ref.@). Therefore, in Fig.l we present 
the results of our calculations for the orbital currents at Q = (j^vr, for three different 
regimes of co-existence of superconductivity and s + id-CDW phase. The gap equation yields 
a maximum of the critical mean field temperature of id-component of the pseudo-gap (TJ) 
(or maximum of entropy) at this Q (see Ref.0). Our numerical calculations show that the 
resulting temperature behavior of the induced magnetic field which is directly proportional 
to the orbital current is indeed slightly differs from the temperature dependence of the s + id- 
CDW gap. One sees from Fig. 1 that the calculated curve reproduces well the observed 
behavior of INS intensity!'! shown in the inset. The values of the hopping integrals were 
chosen (in meV) ti = 100, t 2 = 15 and t 3 = 12. They reproduce well the topology of the 
Fermi surface for underdoped cuprates. 

We also would like to note the following results of our calculations. At T=T C the jump in 
the current strength is reproduced only for the x and y values representing the non-mean field 
character of the coexistence phase between s + ici-CDW and superconductivity. However, 
such a mean-field reduction of the effective gap at T < T c was clearly demonstrated earlier 
by tunneling spectroscopy^. Therefore weather or not a 'jump' at T < T c exist becomes 
an important issue. For example there is no anomalies at T c in /iSR experiments!, but this 
'jump' is clearly visible according to the neutron scattering data!!. 

Let us also comment on the importance of the small s-component of the extended CDW. 
It was shown previously that relatively small s— component of the extended CDW is required 
for the explanation of monotonic increase of the nuclear quadrupolar resonance(NQR)-Cu(2) 



linewidth at T < T c in YBa2Cu306+Jiil. In this connection it is also logical to switch the 
s— component on in our discussion. This component is real and as one can see from Eqs. 
(P^)-(Il) contribute to the orbital current strength if one takes into account the deviation 
of instability vector Q from (71*, n). The values of the s— component was taken as in paperS 
with critical temperature T* = 100K. As one can see the deviation effect together with 
s-component of the pseudogap can reproduce well the observed a 'jump' at T < T c in 
Refs.ll. 

There is an additional argument in favor of the relevance s- component pseudogap with 
respect to the 'jump' of the effective magnetic moment at T < T c . It is connected with the 
dynamical character of the charge- current state. In fact, the discussed internal fields are 



not static. All experiments test the mean squared field yj < h e ^ > that results due to an 
averaging dependent on how fast those fluctuations are. If the s— component exist the energy 
of the sliding (s + id) — CDW condensate depends on the phase of the order parameter Gk 
at T < T c and hence the sliding motion becomes decelerate. Effectively it can be viewed 
as an increase of the measured mean square internal magnetic fields. The situation can 
become even more complicated because, as it was stressed recently the appearance of the 
s— component of CDW leads to a damping of the quasiparticle regimei. In other words 
one can say that the appearance of the s— components of CDW stimulate the nano-scale 
localization phenomenon. The relation of the disorder to the problem of the local magnetic 
fields in YBa2Cu306 +y from experimental point of view was stressed also very recently^. 

Finally we note that the current network pattern, corresponding to the case of Q = 
was discussed by Varmail. We do not touch this case here because according to the neutron 
scattering data the Q is about (vr,7r). Therefore this case seems to be less actual. Further 
experimental studies are required in order to verify the orientation of the observed magnetic 
moments. According to Ref.i they are aligned along c-axis (this is in agreement with our 
discussed (. + W)-CDW scenario), whereas Sidis et al3 reports the polarization of the 
magnetic moments in the copper-oxygen plane. 

In summary, the observed monotonously increasing of the magnetic moment by neutron 



scatteringp'Q is well reproduced our calculations if one identifies the temperature of T* around 
300 K as a critical temperature of (s + id)-CDW phase formation. This value is correlated 
with those T* that were reported earlier as a pseudogap temperature in these compounds. 
We argue that the reported jump in the intensity at T = T c can be attributed to the 
presence of the relatively small s— component of the pseudogap due to deviation of the 
CDW instability vector Q from (tt, 7r). The latter results from the observed topology of 
the Fermi surface and the corresponding behavior of the dynamical charge susceptibility in 
underdoped cuprates0. 
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FIGURES 
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FIG. 1. Calculated intensity of the orbital currents in units of the magnetic fieldllZl 
Hj n4 = J< )+J hcr +J ^ ) ( wri ere r « 2A) at Q = (j^tt, ]jtt) for three different regimes of the 
co-existence between extended CDW and superconductivity as described in the text. Inset shows 
the experimental data taken from Ref. |5[ 
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